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Analytical approximations for (ip 2 ) of a quantized scalar field in ultrastatic asymptotically flat 
spacetimes are obtained. The field is assumed to be both massive and massless, with an arbitrary 
coupling £ to the scalar curvature, and in a zero or nonzero temperature vacuum state. The expres- 
sion for (ip 2 ) is divided into low- and high-frequency parts. The expansion for the high-frequency 
contribution to this quantity is obtained. This expansion is analogous to the DeWitt-Schwinger 
one. As an example, the low-frequency contribution to (p 2 ) is calculated on the background of the 
small perturbed flat spacetime in a quantum state corresponding to the Minkowski vacuum at the 
asymptotic. The limits of the applicability of these approximations are discussed. 
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O . I. INTRODUCTION 

\q • The investigations of black hole evaporation and particle production by an expanding universe have acted as a 
stimulus for a detailed and systematic investigation of the theory of quantum fields propagating on curved spacetimes. 
The main objects to calculate from quantum field theory in curved spacetime are the quantities (v? 2 ) an d {T^) where 
ip is the quantum field and T>f is the stress-energy tensor operator for ip. The renormalized stress-energy tensor 
(Tj^) is an important quantity for the construction of a self-consistent model of an evaporating black hole, while the 
mean-square field (</? 2 ) plays a role in the study of theories with spontaneous symmetry breaking. The functional 
dependence (T^) on metric g M „ allows us to study the evolution of the background geometry driven by the quantum 
fluctuation of the matter fields propagating on it. This is the so-called backreaction, governed by the semiclassical 
Einstein equations 
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G^SttCC). (1) 



JL . However, the exact results for (p> 2 ) and (T^ 1 ) in four dimensions are not numerous (see, for example, Q). Numerical 
qj' computations of these quantities are as a rule extremely intensive 0, 0, 5 . 

One of the most widely used techniques to obtain information about these quantities is the DeWitt-Schwinger (DS) 
L* ■ expansion It may be used to give the expansions for {ip 2 } and (T^) in terms of powers of the small parameter 

^ : — t « i, (2) 

a ; 

where m is the mass of the quantized field and L is the characteristic scale of change of the background gravitational 
field @. 

The analytical approximations to (ip 2 ) and (Ttf) for the conformally coupled massless fields [S El 0> IS IS EH El 
IT^ | give good results. Nevertheless, there still remains the problem of the extension of this type approximations' 
applicability limits. If the quantum field is massive but the mass of the field does not satisfy the condition J2J the 
analytical approximations to (ip 2 ) and (T{f) are even less numerous [5, lil ITTl fl2t IT3L fl4| . 

In this paper, approximate expressions for (<p 2 ) ren of a quantized scalar field in ultrastatic asymptotically flat 
spacetimes are derived. The field is assumed to be both massless or massive with an arbitrary coupling £ to the scalar 
curvature R, and in a zero or nonzero temperature vacuum state. The expression for (<p 2 ) ren is divided into low- and 
high-frequency parts. The Bunch-Parker approach ES] is used for derivation of high-frequency contributions (HFC) 
to these quantities. As in the case of massless field the quantum state of the field with mass m <C 1/1/ is essentially 
determined by the topology of spacetime and the boundary conditions. In this paper such dependence is determined 
by the low-frequency contribution (LFC). As an example, this contribution is calculated on the background of the 
small perturbed flat spacetime in quantum state corresponding to the Minkowski vacuum at the asymptotic. 
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The paper is organized as follows. In Sec. II the expressions for the Euclidian Green's function of a scalar held with 
arbitrary mass and curvature coupling in a ultrastatic spacetime is divided into low- and high-frequency parts. In Sec. 
Ill the WKB approximation of the high-frequency contribution to (ip 2 ) is derived. The low-frequency contribution 
is derived and the renormalization procedure is described in Sec. IV. The results are summarized in Sec. V. In the 
Appendix the expanding of renormalization counterterm for (ip ) in powers of coordinate difference of the separated 
points is described. The units S = c = G = fcs = l are used throughout the paper. The sign conventions are those of 
Misner, Thorne and Wheeler fl(| . 



II. GREEN'S FUNCTION 



The metric in Euclidean section for an ultrastatic spacetime is given by 

ds 2 = dr 2 + g ab dx a dx b , (3) 

where r = it is the Euclidean time, g ab are arbitrary functions of spatial coordinates X 1 , a; 2 , x 3 . (Latin indices run 
from 1 to 3, Greek indices run from to 3.) 

In this paper, the point-splitting method is employed for the regularization of ultraviolet divergences. When the 
points are separated one can show that 

(<p 2 ) U nren=G E (x>*,X»). (4) 

The Euclidean Green's function G E is a solution of the equation 

(□. - m 2 (R(x«)) G^x^- ^T^faf^ ® 

vs ( 'f) 

where m is the mass of the scalar field, £ is its coupling to the scalar curvature R(x a ), g^(x a ) — det g a b(x 1 , x 2 , x 3 ). 
If the scalar field is at zero temperature then S(t, f) and G E {x^ ,x v ) can be expanded as 

1 r + °° 

S(r,f) = — dtoe^- f \ (6) 

2 "" J~oo 



1 f-r°° 

G R (x*,Sr) = — / dLoe^ T -^G E (uj;x a ,x b ). (7) 

If the field is at temperature T, then the Green's function is periodic in r — f with period 1/T. In this case 5(t, f) 
and G E (x 11 , x u ) have the expansions 

S(T,f)=T ]T exp[iw B (T-f)], (8) 



G E (x fi ,x v )=T eM^n{T-f)]G E (LJ n ;x a ,x b ), (9) 



where ui n — litTn. In both cases G e (uj; x a , x b ) satisfies the equation 
1 d ( I , ow " nh , „ s d ~ 



^J g m{x-)g ab {x c )^G E {uj; x a , x b )^j - (c 2 + m 2 + £R) G E (io; x\ ~x b ) 

6^{x a lX b ) 



(10) 



In the case m ^> 1/L, where L is a characteristic scale of the variation of the background gravitational field, it is 
possible to construct the iterative procedure of the solution of this equation with small parameter 1 / mL 0, Il5| . This 
procedure gives the standard expansion of {<p 2 ) in terms of powers of mL. 



In the case 

1 , X 

m<- (11) 
Li 

a small parameter of iterative procedure does not exist. Nevertheless, let us divide G E into low- and high-frequency 
parts, as it was done in 0] 

G E {x»,x v ) = G™ c {x»,x v )+G™ c {x^x v ), (12) 

where 
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G^ c {x^,x u ) = -J dujcos[Lj(T-f)]G E (u};x a ,x"), (13) 

1 r°° 

Gl FC (x fi ,x u ) = - dujcos[uj(T-r)}G E (uj- 1 x a ,x b ) (14) 
n Jl/Xo 

if the scalar field is at zero temperature and 

"o- l 

Gf c (x^,x v )=TG E (0;x a ,x b )+2Tj2 cos[w„(r — f)]G B (u n ; x a , x b ), (15) 

n=l 

oo 

G^ FC (x fi ,x v ) = 2T £ cosK(r-f)]G B (w n ;a: ,s 6 ) (16) 

n— no 

if the field is at temperature T. Then the expansion of G^ FC (x M , x 17 ) in terms of powers of a small parameter 

£WKB = ^ « 1 ( 17 ) 

can be obtained by analogy with the methods of evaluation of the DeWitt-Schwinger expansion. If the field is at 
temperature T then 

Ao = (18) 

and 

£ ™ = 2^L « 1 (19) 
Below the main points of the Bunch and Parker approach |l5( for obtaining G E FC (x ,i , x v ) are outlined. 

III. HIGH-FREQUENCY CONTRIBUTION TO (tp 2 ) 

Let us introduce Riemann normal coordinates y a in 3D space with origin at the point x a 01 • I n these coordinates 
one has 

9ab(y a ) = Vab - \Racbd tfV + ^(y 3 ), (20) 

g ab (y a ) = v ab + \R ab d y c y d + 0(y 3 ), (21) 
5 ( % Q ) = 1 - ^WV + 0(y 3 ), (22) 
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where the coefficients here and below are evaluated at y a = (i.e. at the point x a ), rj a b denotes the three-dimensional 
Euclidean metric. All indices are raised and lowered with the metric rj ab . Defining G(uj; y a ) by 



G(u;;y a ) = JgW(y)G E (u,;y a ) 



(23) 



and retaining in (|10H only terms with coefficients involving two derivatives of the metric or fewer one finds that 
G(w; y a ) satisfies the equation 



ah 



d 2 G 
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.<SG + '-R a c \y c y d 



d 2 G 



QyaQyb 3 c d U » Qya dy b 



C-i)/.' 



G = -5 (3) (y). 



(24) 



Note that quantities R a bcd, Rab, R evaluated in metric © and those, evaluated in 3D metric g a b coincide. Let us 
present 



G(u; y a ) = G q (lo; y a ) + G x {u; y a ) + G 2 (lo- y a ) + 



(25) 



where Gi(u>; y a ) has a geometrical coefficient involving i derivatives of the metric at point y a = 0. Then these functions 
satisfy the equations 



^a 2 G (^;y a ) 



QyaQyb 



co 2 G (co;y a ) = -S (3) (y), 



(26) 



QyaQyb 



^G 1 ( W ;y a ) = 0, 



(27) 



b d 2 G 2 (u:y*) 1 . d d*G (uj;y°) 

V dy«dy» -"G^y) + ~R cd yy 



Go(cu;y a )=0. 



The function Go(w;y a ) satisfies the condition 



b d d 2 G (u ]y a ) b dG (w;y a ) _ 

RcdVV dy«dy» RbV dy* ' 

since Gq(w; y a ) is the function only of r} ab y a y b ^f|. Therefore Eq. i|28|) may be rewritten 



n 



ab d 2 G 2 {uj;y a ) 2T 1 b dG (w;y a ) 

dy*dy» " w G2 ^>y ) + 3 R b y dy a 



G (Lo;y a ) = 0. 



(28) 



(29) 



(30) 



Let us introduce the local momentum space associated with the point y a = by making the 3-dimensional Fourier 
transformation 



G l (to;y a ) 



(31) 



It is not difficult to see that 



G (co;k a ) = 



k 2 



(32) 



Gi(w;fc a ) = 0, 



(33) 
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where k 2 = j] ab k a k b . In momentum space Eg . 1130(1 gives 



(k +uj )G 2 (uj;k ) - -R b k a — 



- (m 2 +£R) G {uo-k a ) = 0. 



(34) 



Hence 



G 2 {w;k) = 



-m 2 -£R 2 R ab k a k b 



(k 2 +LU 2 ) 2 3(fc 2 +CJ 2 ) 3 

Substituting (J2U , ijgjfl. (J33J, ||2H1) in 123 and integrating leads to 

+ 00 

G(^;y a ) 



dkidk 2 dk 3 

3 — exp(ik a y ) 



(2tt) 3 



1 (to 2 + £R) 2 i? ab fc a fc 6 



P+cj 2 (fc 2 +w 2 ) 2 3(fc 2 +w 2 ) 3 



(35) 



1 



87rexp(|w|y) 



-to 2 - (g - l/6)fl 2 _ iWV 1 
|w| y 6y 



(36) 



where y = y r\ a bV a V h ■ Using the definition of G(w; y a ) (|23|l and expansion (|22|l one finds 



G E (u;y a )=[ l + -R ab y a y b )G(uj;y a ) 



87rexp(|w|y) 



-to 2 - (g - l/6)i? | 2 | iUl/V 
|w| y 6 y 



The necessary condition for the validity of this approximation is 

1 1 

W> ^ >> r 

Hence one can evaluate G| FC if the field is at zero temperature 



(37) 



(38) 



1 

7T 
I 

''8^ 



G^ c (Ar,y a )=- / dtu cos[luAt}G e (uj; y a ) 

'l/Ao 



(y 2 + At 2 ) 



\2 
A 



Raby a y 



y\o A 2 (y 2 + Ar 2 ) 6(zj 2 + At 2 ) 

where C is Euler's constant, Ar = r — f . 

If the field is at temperature T then the necessary condition for the validity of expression l|37|l is 



(39) 



n > n » 



2nTL 



and 



Gr c (Ar,y a )=2T ]T cos[w„Ar]G B (w n ; z , x b ) 



(40) 



1 



C+-ln|(y 2 + AT 2 )(27rT) 2 |+V(no) 



~2nT [ n - ~) 2 + (2ttT) 2 (^i 2 - n + ~ 



R ab y a y b 



{y 2 + At 2 ) 6(y 2 + Ar 2 ) J ' 



(41) 



where VK* 2 ) is the logarithmic derivative of the gamma function (i.e., the digamma function). The Plana sum formula 
|18| has been used to compute the sums over n in the last expression. 
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IV. LOW-FREQUENCY CONTRIBUTION TO (ip 2 ) AND RENORMALIZATION PROCEDURE 

The behavior of low-frequency modes is determined by the boundary conditions and the topological structure of the 
spacetime. As an example, let us consider the evaluation of the low-frequency contribution to (y 2 ) on the background 
of the small perturbed flat spacetime in a quantum state corresponding to Minkowski vacuum at the asymptotic. If 
the characteristic scale of the gravitational field inhomogeneity A satisfies the condition 



A 

Y <<:1 



(or XTn Q < 1) , 



(42) 



the low-frequency contributions to ((f 2 ) can be expanded in terms of powers of this small parameter. Below the 
zeroth-order term of this expansion will be used for approximation of the low-frequency contributions to ((p 2 ). This 
means that we choose the long-wave modes approximately coincident with long-wave modes of Minkowski vacuum. 
For these modes in a zero temperature quantum state 



G^ c (Ar,y a ) 



1 



lAo 



(2^) 4 J-x,x 

l/Ao 



, 3 exp (iwAt + ip a y a ) 



(uj 2 + p\ + p\ 



pi 



i 

47T 3 
1 



due iuAT 



l/Ao 

i/A 

-l/Ao 



dp 



P sin(p y) 



y (lu 2 + p 2 + to 2 ) 



due 



exp (~yy/oj 2 + to 2 ) 



In the limit At — > 0, y — > 



1 /-l/Ao 

sin(Ar/A ) 



du>e iuAT 



y 



0(y) 



1 



— to In 




0(y) 



If we take into account conditions and (|17|l . i.e. 



then 



to < — , 



1 + 



TO 
~2~ 



o 

4 



\2 
A 



AgTO 2 



\2 
A 



O (A^to 



The analogous calculations for the temperature quantum state give 

Gl FC (AT,y a )=TG E (0;y a )) + 2T V cos[w n Ar]G B (w n ; y a ) 



n=l 



1 

: 87T 2 " 
TO^ 
2~ 



27rT rip 



[1 - O (T 2 n 2 Ar 2 )] - (2itT) 2 (t 



n ) 



1 - — + In 


4{2irT) 2 n 2 

9 


n 


7U 



+ O (m 2 T 2 nl) 



+ O (\ 2 T 4 n 4 ) 



(43) 



(44) 



(45) 



(46) 



(47) 



The renormalization of {ip 2 ) U nren is achieved by subtracting the renormalization counterterm |l9| and then letting 



U 2 ) = lim [(^ 2 ) 



(48) 
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where 



/DS 



1 1 

: 8tt 2 ct + 8^ 

.2 ! 



5"' 



16?r 2 96tt 2 



(49) 



cr is one-half the square of the distance between the points x and x along the shortest geodesic connecting them, <t m 
is the covariant derivative of cr, and the constant m DS is equal to the mass m of the field for a massive scalar field. For 
a massless field m DS is an arbitrary parameter due to the infrared cutoff in (^ip 2 ^^. A particular choice of the value of 
m DS corresponds to a finite renormalization of the coefficients of terms in the gravitational Lagrangian and must be 
fixed by experiment or observation. The details of the calculations of (4> 2 } DS are discussed in Appendix: 



7r 2 \ y 2 + A r 2 

m 2 Raby a y b 

2 6(y 2 + At 2 ) 



At 2 



(50) 



V. RESULTS 

Using Eqs. ®, (E), JUJ, © and one finds 



i? 



If the field is at temperature T then 

)re«- 167r 2 ^~ 

(2ttT) 2 m ; 



In 



A 2 m 2 s 



L 2 



"-5 



4(2tiT) 2 



+ 2^ (n ) 



48tt 2 8tt 2 



i>{n )- ln(n ) 



2n 



I/ 2 



Let us cite the conditions of the validity of expressions l|51f) , l|52|) once more 



™< z , 



A«A «L or A« 



1 



2nTno 



If also 



no > 1, 



the digamma function ip(no) is given by 



ip(no) = ln(n ) - — + O ( — 
2n \ n o 



Hence the expression (|52[) can be rewritten as 

R ( . I 



!i " 1(>;r 2 P 6 ' ' 1J 



7,2 



48tt 2 



L 2 



+ O (A 2 T 4 ?t,q) 



O (A 2 T 4 n 4 ) 



4(2^T) 2 n 2 , (2ttT) 2 , Q ^ I 



n L 2 



(51) 



(52) 

(53) 
(54) 

(55) 

(56) 



(57) 
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The presence of the arbitrary parameter Ao in the expressions (|51|l is a generic feature of local approximation 
schemes 0, 0, ITol ITTl fl2l IbH l20j . For a conformally invariant field this parameter can be absorbed into the definition 
of constant m DS . 

Note, that the approximation which corresponds to the analytical approximation of Anderson, Hiscock, and Samuel 
for (</? 2 ) in case of ultrastatic asymptotically flat spacetime can be obtained by one use of the high-frequency 
approximation of G B (w; y a ) (see Eq. I|37|l ') for all the values of u>. However the use of the high-frequency approximation 
of G E (oj;y a ) for u> <C 1/Ao does not seem obvious. Nevertheless, in case of a conformally coupled massless field such 
procedure gives good results in the asymptotically flat region. 0, 0, Htl flij . 



Appendix 



In this Appendix the expanding of (</> 2 ) D s m powers of x° 
Let points P(x) and P{x) be connected by the shortest geodesic x a = x a (s), where s is the canonical parameter. 
The functions x a (s) can be expanded in the Taylor series about point P(x) 



1 dxQ As + i^(A s ) 2 -|-i^(A s ) 3 + 0((A s ) 4 ) 



1! ds ' 2! ds 2 v ; 3! ds 3 
where the coefficients are evaluated at P(x). Using the geodesic equation 

d 2 x a + ra dxP dx~< _ Q 
ds 2 ^ 7 ds ds 



one finds 



d 3 x a 
ds 3 



where d 1 denotes the partial derivative with respect to x 1 . Hence Eq. l|Alj) can be rewritten 

x a =~x a + ^u a As + 1 (-Tiyu>) (As) 2 

+^ (-9^0 + 2T» tf r^) u^u^(As) 3 + O ((As) 4 ) , 
where u a = dx a / ds. This equation can be inverted 

u a As = e a + ir>V + (IKsT* + Id^g) 6<W + 0(e 4 ), 



where e Q = x a — x a . If we use determinations ct m and a |2l| , we can write 

a(x,x)=~g af3 u<* U e(As) 2 = \g a pe a ^ + ^ a/3 I>W + ^g^^eW 
+ \g a p (d s T^ 5 ) e?e*e s e e + \g a0 T^ 5 T^ e s e £ e» + O (e 5 ) . 



Hence the resulting expression for (f 2 ) DS is 



'An 2 



e p e'e e e v 



(g a ^y 



m'+[£-~)R 



8tt 2 

i 2 1 R a ae a £ f3 



16tt 2 48tt 2 g a pe a e 



+ 0(6). 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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In coordinates r, y a , where y a are the Riemann normal coordinates with origin at the point x a in 3D space. 

6° = At, e a = y a . (A8) 

At point P in these coordinates also 0] 

r^ 7 - (A9) 

and 

g a p (d v T^ s ) eWe* = (R a cde + R a dce ) y b y c y d y e = 0. (A10) 
After the substitution of these expressions into [|A7J| one finds that (<^ 2 ) DS has the form fBT))) . 
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